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1. Introduction 



The experimental data from HERA on the deep-inelastic scattering (DIS) structure function 
(SF) F 2 [§|-[0], its derivative dlnF 2 /dln(l/;c) g|-[§| and the heavy quark parts F{ c and F% b 



[11] enable us to enter into a very interesting kinematical range for testing the theoretical ideas 
on the behavior of quarks and gluons carrying a very low fraction of momentum of the proton, 
the so-called small-x region. In this limit one expects that the conventional treatment based on the 



Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations []12j does not account for contri- 
butions to the cross section which are leading in a v ln(l/x) and, moreover, the partem distribution 
function (PDFs), in particular the gluon ones, are becoming large and need to develop a high den- 
sity formulation of QCD. 

However, the reasonable agreement between HERA data and the next-to-leading-order (NLO) 
approximation of perturbative QCD has been observed for Q 2 > 2 GeV 2 (see reviews in [|13|] and 
references therein) and, thus, perturbative QCD could describe the evolution of F 2 and its deriva- 
tives up to very low Q 2 values, traditionally explained by soft processes. 

The standard program to study the x behavior of quarks and gluons is carried out by compar- 



ison of data with the numerical solution of the DGLAP equation [12] by fitting the parameters of 
the PDF x-profile at some initial <2q and the QCD energy scale A [?]-[18]. However, for analyzing 



exclusively the low-x region, there is the alternative of doing a simpler analysis by using some of 
the existing analytical solutions of DGLAP evolution in the low-x limit [jT^]-[22]. This was done 



so in [19] where it was pointed out that the HERA small-x data can be interpreted in terms of the 
so-called doubled asymptotic scaling (DAS) phenomenon related to the asymptotic behavior of the 



DGLAP evolution discovered many years ago Q23Q. 

The study of [19] was extended in [2C, 2j_, ^2] to include the finite parts of anomalous dimen- 
sions of Wilson operators 1 . This has led to predictions [2j_ ^2|] of the small-x asymptotic PDF 
form in the framework of the DGLAP dynamics starting at some <2q with the flat function 



fa(Ql) = A a (hereafter a = q,g), 



(1.1) 



where f a are the parton distributions multiplied by x and A a are unknown parameters to be deter- 
mined from the data. 

We refer to the approach of ^"fl , 22] as generalized DAS approximation. In that approach 
the flat initial conditions in Eq. ( |1.1[ ) determine the basic role of the singular parts of anomalous 
dimensions, as in the standard DAS case, while the contribution from finite parts of anomalous 
dimensions and from Wilson coefficients can be considered as corrections which are, however, 



important for better agreement with experimental data. In the present paper, similary to [|19J-[^2j], 
we neglect the contribution from the non-singlet quark component. 

The use of the flat initial condition given in Eq. ( [Ll| ) is supported by the actual experimental 
situation: low-<2 2 data [|24], [I[ f25|, Q] are well described for Q 2 < 0.4 GeV 2 by Regge theory with 
Pomeron intercept Gtp(0) = Ap + 1 = 1.08, closed to the standard (ap(0) = 1) one. The small rise 



'At small x there is another approach based on the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation Jl^], whose 
application will be dicussed below in Appendix A. 

2 In the standard DAS approximation [^] only the singular parts of the anomalous dimensions were used. 
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of HERA data [jl|, ^j, ^] at low Q 2 can be explained, for example, by contributions of higher 



twist operators (see [|22|]). 

The purpose of this paper is to demostrate a good agreement between the predictions from 
the generalized DAS approach and the HERA experimental data [jj], ||] and [0]-[jn|] for SF F2 and 
Fj c and also to compare the predictions for the slope d lnf^/^ \n(l/x) with the HI and ZEUS data 
[@> @> ( see Figs. 1-8). Looking at the HI data points ^ shown in Figs. 5, 6 and 7 one can 
conclude that A(<2 2 ) is independent of x within the experimental uncertainties for fixed Q 2 in the 
range x < 0.01. Indeed, the data are well described by the power behavior 

F 2 (x,Q 2 ) = Cx- X ^\ (1.2) 

where X{Q 2 ) =a\n(Q 2 /K 2 ) with C « 0.1 8, a ss 0.048 andA = 292 MeV %. The linear rise of the 
exponent X{Q 2 ) with ln<2 2 is also explicitly shown in Figs. 5, 6 and 7 by the dashed line. 

The rise of X (Q 2 ) linearly with In Q 2 can be tracted in strong nonperturbative way (see [ 27 ] and 



references therein), i.e., X(Q 2 ) ~ \/a s (Q 2 ). The previous analysis [28], however, demonstrated 
that the rise can be explained naturally in the framework of perturbative QCD. 

The ZEUS and HI Collaborations have also presented [Q, ||] new preliminary data for X(Q 2 ) 
at quite low values of Q 2 . As it is possible to see in Fig. 8 of [Q], the ZEUS value for X(Q 2 ) 
is consistent with a constant ~ 0.1 at Q 2 < 0.6 GeV 2 , as it is expected under the assumption of 
single soft Pomeron exchange within the framework of Regge phenomenology. These points lie 
slightly below the corresponding ZEUS data but all the results are in agreement within modern 
experimental errors. 



It is important to extend the analysis of [21 , 22, 28] to low Q range with the help of well 



known infrared modifications of the strong coupling constant. Indeed, in Ref. [£^], we have used 
the "frozen" and analytic versions (see, [ ]35| ] and [f36|], respectively). 

This contribution is organized as follows. Sections 2 and 3 contain basic formulae, which are 



needed for the present study and were previously obtained in []21|, |22|, |28|, |29|, PCQ . In Sections 4 
and 5 we compare our calculations with HI and ZEUS experimental data and present the obtained 
results. Some discussions can be found in the conclusions. Some preliminary results accounting 
for BFKL corrections in our analysis can be found in Appendix A. It is hoped that the inclusion 
of these corrections will improve the agreement with the experimantal data for F2 and its slope at 
<2 2 ~ 1^2Gev 2 . 

2. Generalized DAS approach 



The flat initial condition ( |l.l| ) corresponds to the case when partem density tend to some con- 
stant value at x — > and at some initial value Qq. The main ingredients of the results ^2| ], 



are: 



Both, the gluon and quark singlet densities are presented in terms of two components (" + " 
and " — ") which are obtained from the analytic <2 2 -dependent expressions of the correspond- 
ing (" + " and " - ") PDF moments. 3 



3 Such an approach has been developed | ^T| ] recently also for the fragmentation function, whose first moments (ie 
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The twist-two part of the " — " component is constant at small x at any values of Q 2 , whereas 
the one of the " + " component grows at Q 2 > Qq as 



~ e 



d++ + \d+ -5- I p 



ln1 ^)' P = 2ln^) 



where a and p are the generalized Ball-Forte variables, 

'"MY 



In 



a s (Q 2 



P = a s (Q ) - a s (Q ), d+ 



12 



412 

27pV 



(2.1) 



(2.2) 



Hereafter we use the notation a s = a s / (47r). The first two coefficients of the QCD /3-function in the 
MS-scheme are j3 = 1 1 - (2/3)/ and j8i = 102 - (1 14/9)/ with / is being the number of active 
quark flavors. 

Note here that the perturbative coupling constant a s (Q 2 ) is different at the leading-order (LO) 
and NLO approximations. Indeed, from the renormalization group equation we can obtain the 
following equations for the coupling constant 

-.2 



1 



a\°{Q 2 ) 



A) In 



A 2 

yv LO 



at the LO approximation and 



J_ ft, 

a s {Q 2 ) + A) n 



A) In 



9L 

A 2 



(2.3) 



(2.4) 



at the NLO approximation. Usually at the NLO level MS-scheme is used, so we apply A = A^- 
below. 



2.1 Parton distributions and the structure function F2 

The results for parton densities and F 2 are following: 
• The structure function F2 has the form: 

F 2 .lo( x ,Q 2 ) = eU,Lo(x,Q 2 ), 
at the LO approximation, where 



(!«?)// 



(2.5) 



(2.6) 



is the average charge square, and 



F 2 (x, Q 2 )=e[ f q (x, Q 2 ) + ^fa s (Q 2 )f g (x, Q 2 



2 

fa(x,Q 2 ) = f+(x,Q 2 )+f-(x,Q 2 ) 
at the NLO approximation. 



(2.7) 



mean multiplicities of quarks and gluons) were analyzed [|32[] . The results are in good agreement with the experimental 
data (see contribution [B3n by Paolo Bolzoni to this Proceedings). 
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The small-x asymptotic results for the LO parton densities f a LO are 



/g,LO ( X 


Q 2 ) 


= (^ + ^/ (a L o)e-^' VL ° 


+ O(plo), 




Q 2 ) 


= i( A s + l A «) P Lo!l ( aL °) e 




fg,ho( x 


Q 2 ) 


= - A -A q e- d - s ^ + 0(x), 




fq,Lo( X 


Q 2 ) 


= A q e- d - SL0 + 0(x), 





where 



d+ = l + 20//(27j3 ), tL = 16// (27 A, 



(2.8) 

(2-9) 

(2.10) 
(2.11) 

(2.12) 



are the regular parts of the anomalous dimensions d + (n) and d- (n), respectively, in the limit 
n — > l 4 . Here n is the variable in Mellin space. The functions I v (v = 0, 1) are related to the 
modified Bessel function I v and to the Bessel function J v by: 



/ v (<7) 



/ v (ff), iis>0 
r v J v (io), i 2 = -l, if s<0 



(2.13) 



At the LO, the variables (7lo and Plo are given by Eq. ( pT| ) when p = 0, i.e. 



Olo = 2WU+ULoln 



Plo 



21n(l/jc) ' 

and the variable slo is given by Eq. (2.2) with fl^°(<2 2 ) as in Eq. (2.3). 
The small-x asymptotic results for the NLO parton densities f£ are 

ft^Q 2 ) 

f g (x,Q 2 ) 
fa(*,Q 2 ) 



(2.14) 



where 



l-d q + _a s (Q 2 )jph{a) + 20a s (Q 2 )I (a) 

A;(Q 2 ,Ql)e- d - Ws - D - p + 0(x), 
A - e -d.{\)s-D. P + 



-d+(\)s-D + p 



D± = d±± - ir d± 
Po 



(2.15) 

+ 0(p),(2.16) 

(2.17) 
(2.18) 

(2.19) 



and similar for D + and D + , 



A g (Q 2 ,Q 2 ) =A g -A+(Q 2 ,Q 2 ). 



(2.20) 



4 We denote the singular and regular parts of a given quantity k(n) in the limit n— > 1 by k/(n— 1) and respectively. 
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The coupling constant a s (Q 2 ) is introduced in Eq. (|L4|). The variables d + , d ++ d + and d- 

are diven in Eqs. (2.2) and ( |2.12 ), respectively. The variables d ++ , d and d\_ have the 

form 



1643 2f 

^ ++ = ^l 36C3 + 33C2 "^2- + T 



with £3 and £2 are Eller functions. 5 



68 ,r 13 / 

46 - 

9 s 243 



,r 23 13 / 
4 ^-18 + 243 



. 4- 



13 f 

23 - 126 (2.21) 



2.2 Effective slopes 



Contrary to the approach in Q19|]-[[22|] vaiious groups have been able to fit the available data 
using a hard input at small x: x , X > with different X values at low and high Q 2 (see [37]-[43]). 
Such results are well-known at low Q 2 values [pq]. At large Q 2 values, for the modern HERA data 
it is also not very surprising, because they cannot distinguish between the behavior based on a steep 
input parton parameterization, at quite large Q 2 , and the steep form acquired after the dynamical 
evolution from a flat initial condition at quite low Q 2 values. 



As it has been mentioned above and shown in y21| , |22| ], the behavior of parton densities and F 2 
given in the Bessel-like form by generalized DAS approach can mimic a power law shape over a 
limited region of x and Q 2 

f a {x,Q 2 )~x-^^ and F 2 (x,Q 2 ) ~ x ~^\ 



The effective slopes X% (x, Q ) and X§^(x, Q ) have the form 



Xf(x,Q 2 ) 



X§f(x,Q 2 ) 



Pi 



f+(x,Q 2 ) J,(g) 
f g {*,Q 2 ) 



I 2 (a)(\-20a,(Q 2 )) + 20a,{Q 2 )h(a)/p 
f q (x,Q 2 ) P h(o)(\-20a s (Q 2 )) + 20a s (Q 2 )I (o)/p ' 
^q ff (x,Q 2 )f+(x,Q 2 ) + (2f)/3a s (Q 2 )Xg ff (x,Q 2 )f+(x,Q 2 



f q (x,Q 2 ) + (2f)/3a s (Q 2 )f g (x,Q 2 ) 
where the exact form of parton densities can be found in ^2||. 



(2.22) 



The results ( |2.22| ) (and also ( |2.23| )-( |2.25[ ) below) are given at the NLO approximation. To 
obtain the LO one, it is necessary to cancel the term ~ a s (<2 2 ) and to use Eqs. (|2~8|)-(pT|) for 
parton densities f a (x, Q 2 ). 

The effective slopes X^ n and Apf depend on the magnitudes A a of the initial PDFs and also 
on the chosen input values of <2q and A. To compare with the experimental data it is necessary the 



exact expressions ( |2.22[ ), but for qualitative analysis it is better to use an approximation. 



5 Note that evaluation of the results (^.8|)-( p!2l| ) need the knowledge of the analytic continuation of the anomalous 
dimansions and coefficient functions. The analytic continuation can be found in Refs. |g4[]. It was used also for the fits 
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2.3 Asymptotic form of the effective slopes 

At quite large values of Q 2 , where the "— " component is negligible, the dependence on the 
initial PD disappears, having in this case for the asymptotic behavior the following expressions 6 : 

h(<y) _ . i 



Xf^(x,Q 2 ) 



Xf^(x,Q 2 ) 



1 eff,as / ^2 



Io(o)~ H 41n(l/x)' 

I 2 (o){l-20a s {Q 2 ))+20a s (Q 2 )h(o)/ P 

h(o)(l-20a s (Q2))+20a s (Q2)I Q (o)/p 



(2.23) 



+ 



\0a s (Q 2 



41n(l/x) pln(l/x)' 



(2.24) 



(2.25) 



r 41n(l/jc) pln(l/x) q < ' pln(l/jc)' 

where the symbol m marks the approximation obtained in the expansion of the usual and modified 



Bessel functions in (|2J_3j). These approximations are accurate only at very large a values (i.e. at 
very large Q 2 and/or very small x). 



As one can see from Eqs. (123) and (124), the gluon effective slope A5 ff is larger than the 



quark slope A| ff , which is in excellent agreement with MRS [i4] and GRV [45] analyses. 

We would like to note that at the NLO approximation the slope Apf' as (x, <2 2 ) lies between 
quark and gluon ones but closely to quark slope A| ff ' as (x, Q 2 ). Indeed, 



Af as (*,e 2 ) - Ap e f' as (*,e 2 ) 



i 



1 eff.as / n 

% 2 \x,Q , 



eff.as 



x,Q 2 ) 



41n(l/x) 

3a s (Q 2 ) 



+ 26a s {Q 2 



1 



2p ln(l/x) 



pln(lA) 



(2.26) 
(2.27) 



Both slopes A f f ff (x, Q 2 ) decrease with increasing x (see Fig. 5). A x-dependence of the slope 
should not appear for PDFs within a Regge type asymptotic (x~ A ) and precise measurement of the 
slope A„ ff (x, <2 2 ) may lead to the possibility to verify the type of small -x asymptotics of partem 
distributions. 



3. F 2 CC and structure functions 



Recently the HI [Q, |] and ZEUS % |lg] Collaborations at HERA presented new data 7 . More- 
over, the preliminary combine HI and ZEUS data of F 2 cc (x, Q 2 ) and F 2 w, (x, <2 2 ) has been demon- 
strated recently (see [U|). 



6 The asymptotic formulae given in Eqs. ( 2.23 )-( 2.25) work quite well at any Q 2 > values, because at Q 2 = Qq 
the values of A„ ff and X^f are equal zero. The use of approximations in Eqs. ( 2.23 )-( 2.25) instead of the exact results 
given in Eq. (2.22) underestimates (overestimates) only slightly the gluon (quark) slope at Q 2 > Qq. 

7 The papers contain also the references on the previous data on deep-inelastic (DIS) structure functions 

(SFs) Fj' and Fy at small x values. 



7 



Small x behavior of partem distributions 



Anatoly Kotikov 



In the framework of DGLAP dynamics [12], there are two basic methods to study heavy- 
flavour physics. One of them [46] is based on the massless PDF evolution of partem and the other 
one on the photon-gluon fusion (PGF) process [47]. There are also some interpolating schemes 
(see Ref. [48] and references cited therein). 

Here we present the results of Ref. [49] were we applied compact low-x approximation formu- 
lae for the SFs F"(x, Q 2 ), with hereafter i = c,b, observed [ |3Q| ] in the framework of PGF process 
at the first two orders of perturbation theory to these new HERA experimental data J7|]-[jl~ij]. We 
show a good agreement between experimental data and the approach which found without addi- 
tional free parameters. All PDF parameters have been fitted earlier [^2|, 29] from F 2 (x, Q 2 ) HERA 
experimental data. 

In the framework of the generalized DAS approach, the SFs F 2 c (x,Q 2 ) and F 2 (x,Q 2 ) have 
the following form 

Fi(x,(f)^M l 2 Jl,(f,n 2 )f g (x,n 2 ), (i = c,b) (3.1) 

where M' 2 g (1 , Q 2 , jJ. 2 ) is the first Mellin moment of the so-called gluon coefficient function C' 2g (x, Q 2 ,jJ. 2 ). 
Through NLO, M 2 , g {l, 2 2 ,£t 2 ) exhibits the structure 



where 



(2), 



21 



|, ^ = Q 2 + 4m 2 .. 



+ ff(at), (3.2) 



(3.3) 



3.1 LO results 



The LO coefficient function of PGF can be obtained from the QED case [59] by adjusting 
coupling constants and colour factors, and they read [|(| ^l]] 



where 



2x{[l - 4jc(2 - c)(1 - jc)]j8 - [1 - 2x(l - 2c) + 2x 2 (l - 6c - 4c 2 )]L(/3)}, (3.4) 

1+P 



4cx 



l-x' 



l-j3 



Performing the Mellin transformation 

M 2 Jn,c) 



f —C 2 , g (x,c) 
Jo x 



(3.5) 



(3.6) 



we find at n = 1 (see [30]) 



with 



M ( 2 ]{hc) 



J(c) = -VbXnt, 



-[1 +2(1 -c)J(c)} 



l-y/b 
1 + Vb' 



1 



l+4c 



(3.7) 



(3.8) 



8 Note that similar formulas work well for (see [g2|) for high-energy neutrino-nucleo scattering where the effective 
value of the Bjorken variable x is very small. 
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3.2 NLO results 

The NLO coefficient functions of PGF are rather lengthy and not published in print; they are 
only available as computer codes [63]. For the purpose of this letter, it is sufficient to work in the 
high-energy regime, defined by x <C 1, where they assume the compact form j64| ] 



<^(x,c)=/3*W(l,c), 



(3.9) 



with 



= £C A [5 + (13-10c)/(c) + 6(l-c)/(c)], *^(l,c) 



(2)/ 



■AC A Mfl(\,c), (3.10) 



where Ca=N for the colour gauge group SU(N), J[c) is defined by Eq. ( piS| ), and 

1 



1(c) = -Vb 



C(2) + | In 2 f - ln(fec) In? + 2Li 2 (-f) 



(3.11) 



where ? is given in (3JS) and Li2(x) = — J Q (dy/y)\n{\ — xy) is the dilogarithmic function. 

As already mentioned above (see the end of Section 2), the Mellin transforms of C^ g (x,c) 
exhibit singularities in the limit 8± — > 0, which lead to modifications in Eq. (3T). As was shown 
in Refs. [4(], 21, f22| ], the terms involving l/8± correspond to singularities of the Mellin moments 
Atf g (n) at n — > 1 and depend on the exact form of the subasymptotic low-* behaviour encoded in 
ff(x,n 2 ). The modification is simple: 



1 

Si 



1 



1 



1 



8 ± f^(x^ 2 )J x y 



dy 



(3.12) 



where x. = x/b. In the generalized DAS regime, the + and — components of the gluon PDF exhibit 
the low-* behaviour (|T|)-(|l2l|). We thus have [fl], ||] 



1 h(a(x)) 



p(x)I (a(x)y 



1 1 

^«ln T , 
S x 



(3.13) 



where o and p are given in (|2. ip. 

Because the ratio f g (x,Q 2 )/f+(x,Q 2 ) is rather small at the Q 2 values considered, Eq. (3T) is 
modified to become 

F><(x,Q 2 ) nM 2 , g (l,n 2 ,Ci)f g (x,H 2 ), (3.14) 

where M2,g(l, At 2 ) is obtained from M2. g (n, i u 2 ) by taking the limit n — » 1 and replacing \/(n — 1) — > 
1/5-f. Consequently, one needs to substitute 



M { M,c)^M { M,c) ( 7 = 1,2) 



(3.15) 



in the NLO part of Eq. (3.2). Using the identity 



l dy 



k(a(x)) Jx y \y 



\n(bc) 



7(c) 



(3.16) 
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Figure 1: F%{x, Q 2 ) as a function of x for different g 2 bins. The experimental points are from HI (open 
points) and ZEUS [||] (solid points) at Q 2 > 1.5 GeV 2 . The solid curve represents the NLO fit. The dashed 
curve (hardly distinguishable from the solid one) represents the LO fit. 



we find the Mellin transform of Eq. (3.9) to be 



M^(l,c) 



R^ihc) (j = l,2), 



(3.17) 



with (I, a) (j = 1 , 2) are given in ( |3.10[ ). The rise of the NLO terms as x — > is in agreement 
with earlier investigations 



4. Comparison with experimental data for SF F2 and the slope Ap 2 



9 Note , that <5 + determines the behavior of the slope of gluon density (see ( 2.22 )) and also mostly the slope of SF Ft. 
The form ( 3.13 ) of 8 + is in full agreement with the results ( 2.22 ) for the asymptotic form of the effective slope of gluon 
density. 
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Figure 2: F2(x 7 Q 2 ) as a function of x for different Q 2 bins. The experimental points are same as on Fig. 1. 
The solid curve represents the NLO fit. The dash-dotted curve represents the BFKL-motivated estimation 
for higher-twist corrections to i*2(x, Q 2 ) (see [^]). The dashed curve is obtained from the fits at the NLO, 
when the renormalon contributions of higher-twist terms have been incorporated. 



Using the results of previous section we have analyzed HERA data for F 2 and the slope 
d \nF 2 /d ln(l /x) at small x from the HI and ZEUS Collaborations [p-[@]. 

In order to keep the analysis as simple as possible, we fix / = 4 and a s (M|) =0.1166 (i.e., 
A^ 4 ) = 284 MeV) in agreement with the more recent ZEUS results 

As it is possible to see in Fig. 1 (see also 22]), the twist-two approximation is reasonable 
at Q 2 > 2 GeV 2 . At smaller Q 2 , some modification of the approximation should be considered. 
In Ref. [22] we have added the higher twist corrections. For renormalon model of higher twists, 
we have found a good agreement with experimental data at essentially lower Q 2 values: Q 2 > 0.5 
GeV 2 (see Figs. 2 and 3), but we have added 4 additional parameters: amplitudes of twist-4 and 
twist-6 corrections to quark and gluon densities. 

Moreover, the results of fits in [O] have an important property: they are very similar in LO and 
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Figure 3: F2(x,Q 2 ) as a function of x for different Q 2 bins. The experimental points are from HI [[!]] (open 
points) and ZEUS \^ (solid points) at Q 2 > 0.5 GeV 2 . The solid curve represents the NLO fit. The dashed 
curve is from the fits at the NLO with the renormalon contributions of higher-twist terms incorporated. The 
dash-dotted curve (hardly distinguishable from the dashed one) represents the LO fit with the renormalon 
contributions of higher-twist terms incorporated. 



NLO approximations of perturbation theory. The similarity is related to the fact that the small-* 
asymptotics of the NLO corrections are usually large and negative (see, for example, a s . -corrections 
[g(J 0] to BFKL kernel 0] 10 ). Then, the LO form ~ a s (Q 2 ) for some observable and the NLO 
one ~ (X S {Q 2 ){\ — Ka s {Q 2 )) with a large value of K are similar, because A 3> Alo 11 and, thus, 
a s (Q 2 ) at LO is considerably smaller then a s {Q 2 ) at NLO for HERA Q 2 values. 

In other words, performing some resummation procedure (such as Grunberg's effective-charge 
method [52]), one can see that the results up to NLO approximation may be represented as ~ 



°It seems that it is a property of any processes in which gluons, but not quarks play a basic role. 
11 The equality of a s (M^) at LO and NLO approximations, where Mz is the Z-boson mass, relates A and Alo : 
A< 4 ' = 284 MeV (as in &) corresponds to A L0 = 112 MeV (see 0]). 
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Figure 4: x dependence of Fz(x, Q 2 ) in bins of Q 2 . The experimental data from HI (open points) and ZEUS 
(solid points) are compared with the NLO fits for Q 2 > 0.5 GeV 2 implemented with the canonical (solid 
lines), frozen (dot-dashed lines), and analytic (dashed lines) versions of the strong-coupling constant. For 
comparison, also the results obtained in Ref . [|2| through a fit based on the renormalon model of higher-twist 
terms are shown (dotted lines). 



a v (<2g ff ), where Q 2 S 3> Q 2 . Indeed, from different studies [|3|, 54, 55], it is well known that at 
small -x values the effective argument of the coupling constant is higher then Q 2 . 

Here, to improve the agreement at small Q 2 values without additional parameters, we mod- 
ify the QCD coupling constant. We consider two modifications, which effectively increase the 
argument of the coupling constant at small Q 2 values (in agreement with [ 53[ 54, 55]). 

In one case, which is more phenomenological, we introduce freezing of the coupling constant 
by changing its argument Q 2 — > Q 2 +M 2 , where M p is the p-meson mass (see [p5[]). Thus, in the 
formulae of the Section 2 we should do the following replacement: 



a s {Q 2 )^a fv {Q 2 )=a s {Q 2 +M 2 p ) 



(4.1) 
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Table 1: The result of the LO and NLO fits to HI and ZEUS data for different low Q 2 cuts. In the fits / is 
fixed to 4 flavors. 





A 8 


A q 


Qo [GeV 2 ] 


X 2 /n.o.p. 


Q 2 > 1.5GeV 2 
LO 

LO&an 
LO&fr. 


0.784±.016 
932± 017 
1.022±.018 


0.801±.019 
0.707±.020 
0.650±.020 


0.304±.003 
339± 003 
0.356±.003 


754/609 
632/609 
547/609 


NLO 

NLO&an. 

NLO&fr. 


-0.200±.011 
0.310±.013 
0.180±.012 


0.903±.021 
0.640±.022 
0.780±.022 


0.495±.006 
0.702±.008 
0.661±.007 


798/609 
655/609 
669/609 


Q 1 > 0.5GeV 2 
LO 

LO&an. 
LO&fr. 


0.641±.010 
0.846±.010 
1.127±.011 


0.937±.012 
0.771±.013 
0.534±.015 


0.295±.003 
0.328±.003 
0.358±.003 


1090/662 
803/662 
679/662 


NLO 

NLO&an. 

NLO&fr. 


-0.192±.006 
0.281±.008 
0.205±.007 


1.087±.012 
0.634±.016 
0.650±.016 


0.478±.006 
0.680±.007 
0.589±.006 


1229/662 
633/662 
670/662 



The second possibility incorporates the Shirkov-Solovtsov idea [3g, 56, 57] about analyticity 
of the coupling constant that leads to the additional its power dependence. Then, in the formulae 
of the previous section the coupling constant a s (Q 2 ) should be replaced as follows: 

Po Q ~ a lo 

at the LO approximation and 



2 



a an (G 2 ) = « v (G 2 )-^^ + ..., (4.3) 

at the NLO approximation, where the symbol . . . stands for terms which have negligible contribu- 
tions at Q > 1 GeV [|36|] 12 . 

Figure 4 and Table 1 show a strong improvement of the agreement with experimental data for 
7*2 (almost 2 times!). Similar results can be seen also in Figs. 5 and 6 for the experimental data for 
A|f (x, Q 2 ) at x ~ 10~ 3 , which represents an average of the x-values of HERA experimental data. 
The top dashed line represents the aforementioned linear rise of k(Q 2 ) with ln(2 2 ). 

So, Figures 5-7 demonstrate that the theoretical description of the small-<2 2 ZEUS data for 
Xpf(x, Q 2 ) by NLO QCD is significantly improved by implementing the "frozen" and analytic 
coupling constants Cffr(<2 2 ) and a an (Q 2 ), respectively, which in turn lead to very close results (see 
also HI 0]). 

I2 Note that in ^| more accurate, but essentially mor e cumbersome approximations of a an (Q 2 ) have been 
proposed. We limit ourselves by above simple form (|4.2|), ( |4.3| ) and plan to add the other modifications in our future 
investigations. 
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Figure 5: The values of effective slope as a function of Q 2 for x = 10~ 3 . The experimental points are 
from HI ^ ^] (open points) and ZEUS [Q| (solid points). The solid curve represents the NLO fit. The 
dash-dotted and lower dashed curves represent the NLO fits with "frozen" and analytic coupling constants, 
respectively. The top dashed line represents the fit from [|[| . 



Indeed, the fits for F 2 (x, Q 2 ) in []22|] yielded Q 2 « 0.5-0.8 GeV 2 . So, initially we had kf?{x, Q 1 ,) 
0, as suggested by Eq. dl.l| ). The replacements of Eqs. (4J_), ( |4.2| ) and (43) modify the value of 
Xp(x,Qq). For the "frozen" and analytic coupling constants CCf r (Q 2 ) and a an (<2 2 ), the value of 
^Jjf (*>2o) * s nonzero an d the slopes are quite close to the experimental data at Q 2 m 0.5 GeV 2 . 
Nevertheless, for Q < 0.5 GeV 2 , there is still some disagreement with the data, which needs 
additional investigation. Note that at Q 2 > 0.5 GeV 2 our results are even better the results of 
phenomenological models [42, |66|]. 

Figure 7 shows the ^-dependence of the slope A|?'(x, Q 2 ). One observes good agreement 
between the experimental data and the generalized DAS approach for a broad range of small-* 
values. The absence of a variation with x of Alf (x, Q 2 ) at small Q 2 values is related to the small 
values of the variable p there. 

At large Q 2 values, the ^-dependence of A|^(x, Q 2 ) is rather strong. However, it is well known 
that the boundaries and mean values of the experimental x ranges |Q] increase proportionally with 
Q 2 , which is related to the kinematical restrictions in the HERA experiments: x ~ 10~ 4 x Q 2 (see 
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Figure 6: Q 2 dependence of X^(x,Q 2 ) for an average small-* value of x = 10~ 3 . The experimental data 
from HI (open points) and ZEUS (solid points) are compared with the NLO fits for Q 2 > 0.5 GeV 2 imple- 
mented with the canonical (solid line), frozen (dot-dashed line), and analytic (dashed line) versions of the 
strong-coupling constant. The linear rise of X^f (x, Q 2 ) with In Q 2 as described by Eq. (2) is indicated by the 
straight dashed line. For comparison, also the results obtained in the phenomenological models by Capella 
et al. p2|] (dash-dash-dotted line) and by Donnachie and Landshoff [BfJ (dot-dot-dashed line) are shown. 



||, 28] and, for example, Fig. 1 of [Q]). We show only the case with the "frozen" coupling 
constant because at large Q 2 values all results are very similar. 

From Fig. 7, one can see that HERA experimental data are close to Xpf(x, Q 2 ) at x ~ 10~ 4 -f- 
10~ 5 for Q 2 = 4 GeV 2 and at x ~ 10~ 2 for Q 2 = 100 GeV 2 . Indeed, the correlations between x and 
Q 2 in the form x e g = a x 10~ 4 x Q 2 with a = 0.1 and 1 lead to a modification of the Q 2 evolution 
which starts to resemble In Q 2 , rather than In In Q 2 as is standard [p8|]. 



5. Comparison with experimental data for SF F| 



We are now in a position to explore the phenomenological implications of our results for SF 
F 2 CC . As for our input parameters, we choose m c = 1.25 GeV in agreement with Particle Data 
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Figure 7: The values of effective slope Ap* as a function of Q 2 . The experimental points are same as on 
Fig. 4. The dashed line represents the fit from [JjJ. The solid curves represent the NLO fits with "frozen" 
coupling constant at x = 10~ 2 and x = 10~ 5 . 



Group [|67|]. While the LO result Eq. (3/7) is independent of the unphysical mass scale jU, the NLO 
formula (3.2) does depend on it, due to an incomplete compensation of the /I dependence of a s (/j.) 
by the terms proportional to ln(/i 2 /<2 2 ), the residual jj. dependence being formally beyond NLO. 
In order to fix the theoretical uncertainty resulting from this, we put \i 2 = Q 2 + 4ra 2 (see (3.3)), 
which is the standart scale in heavy quark production. 

The PDF parameters /Iq, A q and A g shown in (LI), have been fixed in the fits of F2 experimen- 
tal data (see the previous section). Their values depend on conditions chosen in the fits: the order 
of perturbation theory and the number / of active quarks. 



Below Z?-quark threshold, the scheme with f = 4 has been used Q22[ , |29|] in the fits of F2 data. 



Note, that the F2 structure function contains F 2 CC as a part. In the fits, the NLO gluon density and 
the LO and NLO quark ones contribute to F 2 C , as the part of to F2. Then, now in PGF scattering the 
LO coefficient function corresponds in ra — > limit to the standart NLO Wilson coefficient 
(together with the product of the LO anomalous dimension y qg and Infra 2 . /Q 2 ). It is a general situa- 
tion, i.e. the coefficient funstion of PGF scattering at some order of perturbation theory corresponds 
to the standart DIS Wilson coefficient with the one step higher order. The reason is following: the 
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Figure 8: F^ix,^) evaluated as functions of a: with the LO matrix elements (dashed lines) and with the 
NLO ones and with the factorization/renormalization scale ji 2 = Q 2 +4m 2 (solid lines). The black points and 
red squares correspond to the the combine HI ZEUS preliminary data [ |TT] ] and HI data [0, ||], respectively. 



standart DIS analysis starts with handbag diagram of photon-quark scattering and photon-gluon 
interaction begins at one-loop level. 

Thus, in our F 2 CC analysis in the LO approximation of PGF process we should take f a (x,Q 2 ) 
extracted from fits of F2 data at / 
our / = 4 NLO twist-two fit [[~ 



4 and NLO approximation. In practice, in Q49| ] we have applied 
1 of HI data for F2 with Q 2 cut: Q 2 > 1.5 GeV 2 , which produces 



Ql = 0.523 GeV 2 , A g = 0.060 and A q = 0.844. 

Correspondingly, the NLO approximation of PGF process needs the gluon density exracted 
from fits of F2 data at NNLO approximation, which is not yet known in generalized DAS regime. 
However, we see from the modern global fits [|^], that the difference between NLO and NNLO 
gluon densities is not so large. So, we can apply the NLO form (2.11) of f a (x, Q ) for our NLO 
PGF analysis, too. 

The results for F 2 CC are prsented in Fig. 8. We can see a good agreement between our compact 



13 The difficulty to extend the analysis f>3fi to NNLO level is related with an appearence of the pole ~ l/(n— l) 2 
in the three-loop corrections to the anomalous dimension y gg (see J50I , |68|). The pole ~ l/(n — l) 2 violates the Bessel- 
like solution ( 2. 1 1 ) of DGLAP equation for PDFs at low x values with the flat initial condition (1.1). 
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formulas (|3j]), ( ^2| ), ( |377| ) and ( pT7| ) an the modem experimental data for F 2 cc (x,Q 



2\ 



structure function. To keep place on Fig. 8, we show only the HI [0, [§] data and the combine 
H1ZEUS preliminary [[Tl]] one. 

The good agreement between generalized double-asymptotic scaling DAS approach used here 
and F2 and F 2 data demonstrates an equal importance of the both partem densities (gluon one and 
sea quark one) at low x. It is due to the fact that F2 relates mostly to the sea quark distribution, 
while the F 2 CC relates mostly to the gluon one. Dropping sea quarks in analyse ledas to the different 
gluon densities extracted from F2 of from F 2 (see, for example, J70[]). 

6. Conclusions 



We have shown the Q 2 -dependence of the structure functions Fi and F 2 CC and of the slope 
XfjF = d\nF2/d\n(\/x) at small-x values in the framework of perturbative QCD. Our twist-two 
results are in very good agreement with precise HERA data at Q 2 > 2 GeV 2 , where perturbative 
theory can be applicable. The application of the "frozen" and analytic coupling constants (Xi r (Q 2 ) 
and cc an (Q 2 ) improves the agreement with the recent HERA data [Q, for the slope Xf 1 (x, Q 2 ) 
for small Q 2 values, Q 2 > 0.5 GeV 2 . 

We presented a compact formula for the heavy-flavour contributions to the proton structure 
functions F2 valid through NLO at small values of Bjorken's x variable. Our results agree with 
modern experimental data [f7|]-[jl~ij] well within errors without a free additional parameters. In 
the Q 2 range probed by the HERA data, our NLO predictions agree very well with the LO ones. 
Since we worked in the fixed-flavour-number scheme, our results are bound to break down for 
Q 2 3> 4m?, which manifests itself by appreciable QCD correction factors and scale dependences. 
As is well known, this problem is conveniently solved by adopting the variable-flavour-number 
scheme, which not considered here. 

As a next step of investigations, we plan to perform combined fits of the Hl&ZEUS data [||] 
ofF 2 (x,Q 2 ), the Hl&ZEUS data @ of F 2 cc (x,Q 2 ) and the HERA data of F 2 bb (x,Q 2 ), using 
the "frozen" and analytic coupling constants in both the LO and NLO approximations, in order 
to improve the agreement with HERA data at small Q 2 values. Several versions of the analytical 
coupling constant will be used. 
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7. Appendix A 



Here we give a short introduction of possible accounting for BFKL corrections to our analysis. 
As it was shown in Ref. [ |2~l| ] for the first two orders of the perturbation theory, it is conveninet to 
start with Mellin moment reprentation. 

7.1 Mellin moment form 

In the following we resume the steps we have followed to reach the small x approximate 
solution of DGLAP shown above (see also [21]) 14 : 

• Use the n-space exact solution for "±"-components 



A±exp 



a s {Ql 



with 



"M) y±(n,a s ) 
A) 



'Ate 



-d±(n)s 



d±(n) 



(Al) 



(A2) 



Expand the perturbatively calculated parts (of anomalous dimensions and coefficient func- 
tions) in the vicinity of the point n = 1 . 



14 To work with BFKL formulas in the most symmetric way, in this Appendix we will use the normalization of the 
anomalous dimensions deviated by the factor "— 1 /2" from the DIS standard notation. 
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The singular part of the "+"-component with the form (hereafter n = 1 + ft)) 

A a (o k e- lsLo/w (A3) 

leads to modified Bessel functions in the x-space in the form 

(*+l)/2 



I k+ Al^/ ds\nx\ (A4) 



The regular part B(n)exp (— d{n)s) leads to the additional coefficient (see Ref. [21] and 
Appendix there) 



B{\)e- d ^ s + 0{Jds/\nx) 



behind of the modified Bessel function ( |A4[ ) in the x-space. Because the accuracy is 0(y d + s/ lnx), 
it is necessary to use only the basic term of Eq. ( A4), i.e. all terms co k in front of exp (— d + / ft)), 



with the exception of one with the smaller k value, can be neglected. 

If the singular part at n — > 1 is absent (as in the case of the "—"-component), i.e. Jl = in 



(|A3|), the result in the x-space is determined by B{\)exp(—d(l)s) with accuracy 0(x) 



7.2 BFKL corrections 

We would like to stress that the applicability of the above recipe (to constract the small-x 
solution which was shown in the previous subsection) is not limited by the order in perturbation 
theory but by the form of the singular part of the anomalous dimensions. At the first two orders of 
perturbation theory the singular part is proportional to ~ ft) 1 but this behaviour does not remain 



at higher orders. The most singular terms have been calculated in [14]. For example, the singular 



part of the "+"-component of the anomalous dimension matrix has the following form 

^ k 

CO 



Y+(co,a s ) = y(o), a s .) + O ( a s ( — ) ) (A5) 



where the terms ~ 0[a s {a s / co) ) have been evaluated in Ref. [pOQ. 



The BFKL anomalous dimension /(ft), a) is obtained by solving the implicit equation 

i = ^(tf..*)), 

where the characteristic function %(y) has the following expression in terms of the Euler *F- 
function: 

X (r) = 2w(i)-w(r)-w(i- 7 ), y( y ) = J(ln ™ } 
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7.2.1 Expansions 

The expansion of x(n, /) in powers of 7 gives: 

1 °° 
1 k=\ 

that can be rewritten as the following exact relation for 7 itself 



(A6) 



7=7o 



l+2£C(2^ + 3)7 2 ' 



jfc+3 



k=\ 



7o 



4Caci s 
ft) 



(A7) 



where 70 is the singular part of the LO gluon-gluon anomalous dimesion. Solwing above eqution 
by interations, we have 



7=7o + £c^ 



+1 



k=3 



where the new coefficients are 



C 3 =2C(3), C 4 = 0, C 5 = 2C(5), C 6 = 12C 2 (3), C 7 = 2^(7), ... 



(A8) 



Incorporating the BFKL term to the renormalization exponent leads to the following replace- 



ment 



exp{- 



1 f<Q 2 ) da f+(a) 



1 r a AQ 2 ) da 
}^exp{-— / y(a)} 



PoJa s (Q 2 ) a 2 CO 

The contribution of the additional term L~=3Q:7o +1 i n the r.h.s. of (4.8) has the following 



form 



PO J a k=3 P° a k=3 k 



(A9) 



and, thus, it has additional factor k in the denominator. So, it gives a hopeness that in the form ( |A9j ) 
the BFKL contributions will be not so large as usual. 

7.2.2 Exact contribution 



Now we considere the BFKL contribution to the r.h.s. of ( A9) without any axpansions. Using 



integration by parts procedure, we obtain the following expression 



f da , . 1 . . f da dy(a) 1 . , f dy 

/ — Y (a) = —Y{a)+ / (±-L = — 7 ( a )+ / -L , 

J a a J a da a J a 



(AlO) 



where the inverse coupling constant is proportioanl to the characteristic function %(y) in (A6) 



1 4C A 

- = — X 7 ): 

a CO 
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Thus, the last integral in ( A10) can be evaluated exactly as 

r(l-y) 



2»P(l)y+ln- 



r(y) 



So, the needed contribution in the r.h.s. of (A9) can be represented in the following form 



1 [<Q 2 )da , 
PoJaAQl) a 1 



R(Q 2 ) 
R{Ql) ' 



(All) 



where the R value is (see also Ref. [72]) 



R(Q 2 ) 

with the new parameter do 



r(y) 
r(l-y) 



do 



ex P {-j y(^(y)+^(i-y))} 



(A12) 



do 



4C A 
poco 



a) 



When a s — >• (i.e. in the considered case y — > 0) we recover the singular part of the LO 
contributions (see the previous section) 

R(Q 2 ) -> y^ -> ^ rfo 



The transform of the Mellin moments in the form ( |A11[ ) and ( |A12[ ) to the Bjorken x-space, it 
is not a trivial problem. Author plans to return to this problem in his future work. 
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